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Abstract. The effect of strong magnetic fields parallel to GaAs/GaxAli_xAs 
interface on the subband structure of a 2D electron layer is investigated the- 
oretically. The system with two levels occupied in zero magnetic field is 
considered and the magnetic field induced depletion of the second subband 
is studied. The confining potential and the electron energy dispersion re- 
lations are calculated sclf-consistently, the electron - electron interaction is 
taken into account in the Hartree approximation. 



Recently, there was suggested ^ that the deviations of the 2D Fermi hne 
from the circular shape, which are due to the combined influence of an ap- 
proximately triangular potential well and of the parallel magnetic field, may 
play an important role in the theory of magnetotransport in two-dimensional 
inversion layers at interfaces of GaAs/AlxGai_xAs heterostructures. 

Zawadzki, Klahn and Merkt |§] studied theoretically Fermi lines of narrow- 
band-gap semiconductors of the In-Sb type. Their analysis is based on two 
main simplifications: i) The confining potential Vconfiz) is taken in the form 
of a triangular potential well, ii) The interface is considered as an unpenetra- 
ble hard wall which is accounted for by an appropriate boundary condition 
for the wave function. In this model the spacing of levels due to the mag- 
netic quantization is inversely proportional to the effective mass m, while 
the spacing of subbands due to the confining potential is proportional to 
Thus, for narrow-gap semiconductors with small effective masses 
the electric and magnetic effects become comparable for not to high mag- 
netic fields and for this range of fields also the change of Fermi lines from 
the circular to the 'egg-like' form is expected. 

In semiconductors with larger effective mass like GaAs the effect of par- 
allel magnetic field should be less pronounced and, therefore, it was till 
now treated mainly by the perturbation theory j^. On the other hand it 
is well known that due to the small conduction band offset between GaAs 
and AlxGai_xAs the interface wall is rather soft and electrons can partly 
penetrate into AlxGai_xAs. Thus, to obtain quantitatively better results, 
it is desirable to go beyond the triangular well approximation or the pertur- 
bation approach. For this reasons we decided to performe the full numerical 
self-consistent study of the influence of the in-plane magnetic field on the 
subband shape of 2D electron gas confined to GaAs/AlxGai_xAs interface 
and this contribution presents its results. 

It is now well-established that the numerical self-consistent calculation 
based on the effective mass approximation describes correctly the observed 
electron subband structure in the zero magnetic field. Moreover, generally 
accepted methods of such calculation including the Hartree approximation 
should be valid in the presence of the magnetic field as well as in the zero 
field case. The self-consistency requirement, i.e. the requirement of solving 
coupled Poisson and Schrodinger equations, is an important aspect of the 
energy spectra calculations in doped heterostructures as the charge distri- 
bution reacts on the confining potential which itself determines the charge 
distribution. In strong magnetic fields there is one more reason for the 
self-consistency, the charge redistribution caused by the magnetic field. 



1 



For doped GaAs/AlxGai_xAs heterostructures the total charge density 
q{z) entering the Poisson equation can be spUtted into parts correspond- 
ing to concentrations of electrons, their parent donors in AlxGai_xAs and 
ionized residual acceptors in GaAs: 



Nt{z)+N-{z) 



(1) 



g{z)=e[Ne{z)-..a 

We accept a usual approximation of constant impurity concentrations and 
assume donors and acceptors to be ionized within certain finite intervals 1^ 

and la. N+{z) = for -ld-w< z< -w and N-{z) = TV^ for < z < /„, 
w is the spacer thickness (see figure 1). The confining potential 



Vconfiz) = Vb{z) + Vs.c. (2) 

is a sum of the step function Vb{z) = VbQ{—z) corresponding to the conduc- 
tion band discontinuity and of the Hartree term Vs_c.{z) determined from 
the Poisson equation 



(3) 



dz^ £ 

The conduction band offset Vb and the dielectric constant e enter our calcu- 
lations as input parameters. 

The simplest semiempirical model working quantitatively for the lowest 
conduction states of GaAs/AlxGai_xAs heterostructures is used to solve the 
Schrodinger equation in the envelope function approximation. The envelope 
function is assumed to be built from host quantum states belonging to a 
single parabolic band. The eff'cct of the eff'ective mass mismatch is com- 
pletely neglected and the envelope functions of GaAs and AlxGai_xAs are 
smoothly matched at the interface. 

Due to the translational invariance in the layer plane the wave function 
V'Q.(r) can be factorized 
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and the Schrodinger equation may be written as 



(4) 



2m dz"^ 2m ^ 



e I Byz) - eVconf {z) 



2m 



^i,k^ (^) = 

(^) 



(5) 
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where the magnetic field parallel to the layer plane x — y has a form B = 
(0, By,0). As the heterostructure is in an electric and thermodynamic equi- 
librium two additional conditions must be fulfilled: the charge balance con- 
dition and the constant chemical potential condition. For details about the 
self-consistent procedure for the zero magnetic field (but relevant for B ^ 
as well) we refer to works by Stern and Das Sarma Q and Ando From 
the latter paper we also took values of the band offset Vb = SOOmeV and the 
dielectric constant e = 12.9. 

To distinguish between the charge redistribution resulting from the stan- 
dard self-consistent loop in zero magnetic field and changes due to the mag- 
netic field we performed the self-consistent procedure in two steps. First, 
the electron structure was calculated self-consistently for 5 = and Ne{z) 
and Vconf obtained from this calculation were used as input values for the 
electron structure determined for B ^ 0. In this way the 'intermediate' 
results {B ^ 0, but self-consistency only for B = 0) were obtained. In the 
second step the procedure continued by the full self-consistent calculation 
for B ^ yielding the final resuls. 

The parameters Nd = 2 x lO^^cm^^, Na = lO^^cm"^ and w = 2nm 
were chosen to obtain the electron system of A'^e ~ 14 x 10^^ cm^^ having 
two levels occupied in zero magnetic field and depleting the second level 
at ~ lOT. Note that in this case the 'intermediate' results are all based 
on the electron density corresponding to two occupied levels while the full 
self-consistent study includes the magnetic field induced transfer of electrons 
from the second subband to the first one into the solution of coupled Poisson 
and Schrodinger equations. 

In the zero magnetic field the in-plane electron motion and its out-of- 
plane component along the z-axis are completely independent. It means 
that all electrons within a subband are described by a single localized wave- 
function cpi (z) regardless their energies Ei{kx) and wavevectors k^, ky. In 
our case, with two occupied subbands, we have two different wavefunctions; 
their centres of mass {z)q and {z)i determine the averaged distances of elec- 
trons in subbands from the interface. 

When the magnetic field is applied the effective electro-magnetic poten- 
tial Veff composed from the harmonic magnetic potential and the confining 
potential Vconf is built 

^e// = - ^of - eVconfiz). (6) 

The centre zq of the magnetic part of the effective potential is related to the 
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wave vector component kx by zo = hkx/mu). Thus the magnetic field couples 
the electron motion in x and z directions and for each kx new (^j {z) must 
be calculated. Also the energy spectrum Ei{kx) will deviate from the original 
parabolic dependence on kx and Ei[kx) ^ Ei{ — kx) due to the breakdown 
of the time reversal symmetry. This results in the asymmetric Fermi lines. 
New eigenfunctions fi^k^ (z) are shifted from their original positions {z)o 
and {z)i obtained for B = and, therefore, also the charge distribution 
described by their squares is changed. 

There exists a relation between the centre of mass of the wave function 
V^i,kx (^) shape of the energy spectrum curve Ei(kx) 



which makes it possible to calculate this quantity without numerical diffi- 
culties. As there is one to one correspondence between {z)i^kx ^^'^ kx and 
as we are interested mainly in the charge redistribution induced by the in- 
plane magnetic field. Figure 2 presents the energy subbands as functions of 
instead of kx, together with the shape of the self-consistent potential 
Vconf- The dashed lines correspond to 'intermediate' results while the full 
lines describe the results of the full self-consistent calculations. 

We already mentioned that in the zero magnetic field electrons reach all 
possible energies for just two centres of mass {z)o and {z)i. If the magnetic 
field is applied the electrons with low energies remain locahzed approxi- 
mately around these two values but the electrons with energies closer to the 
Fermi level are shifted either to the interface or deep into the bulk GaAs, 
depending on the direction of their motion. The minimum and maximum 
distances {z)i^kx,min and {z)i^k^^jnax are reached for the Fermi energy and cor- 
respond to electrons moving in opposite directions. Note that the velocity 
is related to the energy spectrum by 

The electrons in the second subband are more sensitive to the magnetic field 
and the second subband is emptied at B = 11. 8T for 'intermediate' and at 
B = 12T for full self-consistent calculation. 

The difference between two types of calculations increases with increas- 
ing magnetic field. This is particularly valid for depressions of the confining 
potential far from the interface which appear in results of full self-consistent 




1 dEi{kx) 



(7) 



fiuj dk. 
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calculations. Note that in spite of this the maximum distance of an elec- 
tron from the interface is not reached for the maximum magnetic field, but 
approximately for S = 12T. 

Figure 3 shows the Fermi lines based both on 'intermediate' and full 
calculations. Two concentric circles correspond to two parabolic subbands 
in the zero magnetic field. With increasing field the area of the second 
subband Fermi surface decreases and the Fermi line of the first subband 
takes the 'egg-like' form. Its left half corresponding to electrons close to the 
interface remains approximately circular, while the right half, describing the 
electrons which move far from the interface in the bulk GaAs takes nearly 
parabolic form. The gross qualitative features of both 'intermediate' and full 
resuls are very similar. The closer look shows marked differences between 
the full and dashed lines. Because of the large magnetic field induced charge 
redistribution, the full lines are qualitatively different having rather the 
'pear-like' than the 'egg - like' shapes. 

Individual wave functions are even more influenced by the magnetic field. 
Figure 4 shows the first subband wave functions (fo^kx i^) corresponding to 
{z)o,kx,min and {z)o,kx,max- In Comparison with the case B = the function 
corresponding to {z)o^kx,min is narrower and shifted to AlxGai_xAs spacer 
and dotted region, while the function corresponding to {z)Q^kx,max is broad 
and almost entirely inside the bulk GaAs. Note a double-peak structure of 
the wave function for B = 6T. 

We conclude that as in the case of the zero magnetic field the electron 
structure of GaAs/AlxGai_xAs heterostructure in strong in-plane fields has 
to be calculated self-consistently if we are interested in semi-quantitative 
results. Both the energy spectra and the confining potential are affected by 
the charge redistribution caused by the magnetic field. On the other hand 
the total number of electrons remain constant. It is due to the fact that 
the highest density of states corresponds to the bottom part of the energy 
spectra which are less influenced by the self-consistency. 
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Figure captions 



Figure 1. The charge distribution in a single modulation-doped GaAs/AlxGai_xAs 
hetero junction with quasi-2DEG at the interface in the depletion length 
model. 

Figure 2. Electron cigcncncrgies E(z) calculated as a function of {z)i^kx 
and the confining potential V{z). The energies are measured from the Fermi 
energy. Dashed lines correspond to the 'intermediate' results, full lines to 
the results of the full self-consistent calculations. 

Figure 3. Lines of constant Fermi energie in (kx, ky) space are shown for the 
two-subband system depleting the second subband at 11. 8T ('intermediate' 
result) resp. 12T (full self-consistent result). The distances in A;-space are 
measured from their centres. 

Figure 4. Wave functions of first subband electrons with (2;)o,fc^,mm) {z)o,kx,max 
illustrate the charge redistribution due to in-plane magnetic fields. 



